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Abstract 

The possibility that both the baryon asymmetry and dark matter arise from the 
late decay of a population of supersymmetric particles is considered. If the decay 
takes place below the LSP freeze out temperature, a nonthermal distribution of 
LSPs results. With conserved R parity these relic LSPs contribute to the dark 
matter density. A net asymmetry can exist in the population of decaying particles 
if it arises from coherent production along a supersymmetric flat direction. The 
asymmetry is transferred to baryons if the condensate decays through the lowest 
order nonrenormalizable operators which couple to R odd combinations of standard 
model particles. This also ensures at least one LSP per decay. The relic baryon 
and LSP number densities are then roughly equal. The ratio of baryon to dark 
matter densities is then naturally f^/f^LSP ~ C(m&/mLSp)- The resulting upper 
limit on the LSP mass is model dependent but in the range (9(30 — 140) GeV. The 
total relic density is related to the order at which the flat direction which gives rise 
to the condensate is lifted. The observed density is obtained for a direction which 
is lifted by a fourth order Planck scale suppressed operator in the superpotential. 
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1. Introduction 

The baryon asymmetry and dark matter density may provide indirect win- 
dows to very early epochs in the evolution of the universe, and to physics at large 
energy scales. In most scenarios the physical mechanisms which give rise to the 
baryon asymmetry and dark matter are unrelated. For example, in supersymmet- 
ric theories the dark matter density is usually assumed to result from the freeze 
out of the lightest supersymmetric particle (LSP). If R parity is unbroken the LSP 
is stable, and the relic LSPs make up the dark matter. The baryon asymmetry 
is usually assumed to arise either at the electroweak phase transition [1], by the 
Affleck-Dine mechanism in which a coherent condensate carrying baryon number 
is generated along a supersymmetric flat direction [2], or the out of equilibrium 
decay of massive particles through baryon and CP violating interactions [3]. In 
all these mechanisms the dark matter and baryon densities are a priori unrelated. 
This is not surprising since the LSP carries a multiplicative quantum number while 
baryon number is additive. The processes which lead to the respective relic densi- 
ties are therefore distinct. That the baryon and dark matter densities are in fact 
the same within a few orders is not necessarily a direct consequence of any of these 
mechanisms, and seems fortuitous. 

Here I suggest an alternate supersymmetric mechanism in which both baryons 
and dark matter arise from the late decay of a weak scale mass particle. As dis- 
cussed below if the mass of the decaying particle is above the LSP mass, and 
the population of decaying particles carries a large asymmetry, then (optimally) 
roughly equal numbers of baryons and LSPs result from the decay. If the temper- 
ature at the era of decay is low enough, the LSPs do not rethermalize, and the 
relic density is determined by the decaying particle density. The ratio of baryon to 
dark matter density in this scheme is then proportional to the ratio of the lightest 
baryon mass to LSP mass, Q^/Q^p ~ Oim^/m^p). For an LSP with weak scale 
mass, this gives roughly the correct ratio, fi^/fiLSP ~ O(10^ 1 — 1CT 2 ). This result 
is reminiscent of the analogous relation in technicolor theories if the lightest tech- 
nibaryon makes up the dark matter. There the electroweak anomaly ensures that 
at high temperatures the baryon and technibaryon number densities are roughly 
equal [4]. Here however, the LSP density is protected from erasure by the low 
temperature at the time of decay, rather than an additive quantum number. 
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In order for this mechanism of relating the baryon and dark matter densities 
to be operative the decay must occur below the LSP thermalization temperature, 
but above the temperature at which nucleosynthesis takes place. This can happen 
if the decaying particle is coupled to standard model fields by nonrenormalizable 
operators suppressed by an intermediate scale, somewhat below the GUT scale [5]. 
These operators must carry baryon number if any asymmetry is to result, and be 
odd under R parity if at least one LSP is to result from each decay. In addition, 
there should be a large particle-antiparticle asymmetry in the decaying population 
if order one baryon per decay is to result. Such large asymmetries can result from 
the coherent production of scalar fields along supersymmetric flat directions. Flat 
directions are likely to be generic features of supersymmetric theories. Finally, in 
order that the total density of the universe have the observed value now, the number 
density of the late decaying particles should be less than thermal at the time of 
decay. Far too many LSPs would remain if the decaying particles had thermal 
number density. A subthermal number density in fact naturally occurs for coherent 
production along flat directions which are lifted by Planck scale suppressed terms 
in the superpotential [6]. The density in the condensate, and therefore the total 
relic density, is related to the order at which the flat direction is lifted. All the 
ingredients for this late decay scenario therefore exist in supersymmetric theories. 

2. Requirements for Baryons and LSPs from Late Decay 

A number of requirements must be met if the late decay scenario for the origin 
of the baryon asymmetry and dark matter is to be realized within supersymmetry. 
In most SUSY models the LSP is typically a neutralino, a linear combination of 
gaugino and Higgsino. If the relic LSPs are to act as dark matter, they must be 
stable as the result of some symmetry. Since the neutralino is Majorana, this must 
be a discrete symmetry, giving a multiplicative quantum number. In what follows 
I will assume the required symmetry is R parity. If the decaying particle is much 
heavier than the LSP then multiple LSPs can in principle be produced in the decay 
chain. However, as discussed below at most one unit of baryon number can result 
from each decay. So unless the LSP is very light, there should not be too many LSPs 
per decay. In order to guarantee that at least one LSP results from each decay, 
the decaying particle should be odd under R parity from the low energy point of 
view. If the mass of the decaying particle is in the range mLSP < Tn^ < 2mLSP, 
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then precisely one LSP results per decay. For simplicity this will be assumed to be 
the case. The decaying particle then also has weak scale mass. 

If the number density of relic LSPs is to be determined by the density of de- 
caying particles, the temperature during the decay epoch should be less than the 
LSP equilibration temperature. If the decay takes place above this temperature, 
the relic LSP density is determined by freeze out, as in the usual scenario. For 
particles with weak scale annihilation cross section and mass, the equilibration 
temperature is roughly T ~ ^j^lsp- With the LSP mass in the range discussed 
below this corresponds to roughly T ~ (9(1 GeV). In addition to this upper limit 
on the temperature at the time of decay there is a lower limit arising from nucle- 
osynthesis. If decays take place during or after nucleosynthesis the light element 
abundances can be modified by photodissociation and photoproduction by decay 
products [7,8]. This can be avoided for T ^ 1 MeV since the weak interactions are 
in equilibrium and the usual neutron to proton ratio results. The decay tempera- 
ture must therefore lie in the window 1 MeV ^ T ^ 1 GeV. The decay rate, T, and 
decay temperature, are related by Tj ~ y^90/g^7r 2 TM p , where M p = m p / a/87t 
is the reduced Planck mass, and g* is the effective number of degrees of freedom 
(g* ~ 10.75 for T ~ 1 MeV). With weak scale mass, such a slow decay rate implies 
the decaying particle must couple to standard model fields only through nonrenor- 
malizable interactions. Decay through Planck scale suppressed couplings leads to a 
decay temperature much too low to avoid the bounds from nucleosynthesis [9,10]. 
However, a decay temperature of order the nucleosynthesis bound in fact results if 
the particle decays through dimension 5 operators suppressed by a scale somewhat 
below the GUT scale [5] . For the 3-body decays discussed below 



where A 2 = |A| 2 is a sum over generations in the final state, X/M is the coefficient 
of the operator, is the mass of the decaying particle, and final state masses have 
been neglected for simplicity. This gives a decay temperature of 



A decay temperature in the window given above can be obtained for 3 x 10 1U GeV 



r ~ 



(8tt) 3 M 2 



(1) 
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^ Mj A ^ 3 x 10 13 GeV. Although this is probably too low to be associated directly 
with the GUT scale, it could arise from an intermediate scale. 

Producing a baryon asymmetry in the decay imposes a number of additional 
requirements. The particle must of course decay through an operator which trans- 
forms under U(1)b with respect to the standard model fields. In principle non- 
renormalizable couplings could arise from D type Kahler potential terms or F type 
superpotential terms. However, with conserved R parity, the gauge invariant op- 
erators which carry baryon number contain at least 3 standard model fields. A 
Kahler potential coupling of this type to the decaying particle is dimension 6, but 
a superpotential coupling to 3 fields is dimension 5. The only invariant made out 
of 3 standard model fields which carries U(1)b is udd. The unique superpotential 
coupling which satisfies the requirements is therefore 

W = A <fmdd (3) 

where <fi is the decaying particle and generation indices are suppressed. Notice that 
udd is odd under R parity. So with an unbroken R parity (at least) one LSP results 
from each decay. In addition, if R parity is to remain unbroken after the decay, 
= must be the ground state. 

Depending on the specific model might decay through other dimension 5 
terms in addition to (3). Decay through superpotential couplings to the other R 
odd combinations of 3 standard model fields, namely LdQ, and LLe (which do 
not carry baryon number) would still give at least one LSP per decay, but dilutes 
the baryon number (for the decay from a condensate with a particle-antiparticle 
asymmetry discussed below). In GUT theories, such operators are in general re- 
lated by GUT symmetries. For example, in 577(5) models udd, LdQ and LLe are 
contained in 5510. The existence of these other decay channels related by SU(5) 
would dilute the baryon number by a factor | . All other dimension 5 couplings are 
through operators which do not carry baryon number. These couplings include: 1) 
superpotential couplings to R even combinations of standard model fields, namely 
QH u u, QH^d, and LH^e, 2) Kahler potential couplings 

M^xh (4) 
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where x is a light field, and 3) dependence of the gauge kinetic functions 

T 0\V"\Vr, (5) 



32tt 2 M 



where W a is the field strength for a light gauge supermultiplet. All of these decay 
modes of course do not contribute to the baryon asymmetry. However, if < 
2mLSP n o LSPs result either. So if is light enough these decay modes do not 
affect the relic f^/^LSP- Finally, the coupling 4>LH U) if present, would allow decay 
through a renormalizable operator, giving a very large decay temperature. In 
addition, it would cause to pair up with some linear combination of neutrinos 
after electroweak symmetry breaking, giving a Dirac neutrino with weak scale mass. 
This (dangerous) coupling must therefore be restricted in some way (as in the toy 
model given in the next section). 

The decay through the operator (3) can in principle lead to a net baryon 
asymmetry, parameterized by e = (iV&)/ (A^lsp) , where (iVj,) and (Alsp) are the 
average number of baryons and LSPs resulting from each decay. In order for 
^&/^lsp ~ CK m &/ m LSp) to hold, e should not be too small. Direct production 
of a baryon asymmetry in the decay requires decay channels which carry different 
baryon number, final state interactions, and CP violating interference terms which 
contain at least two baryon violating couplings [11]. With conserved R parity, 
baryon number is violated only by nonrenormalizable operators, giving negligi- 
ble interference terms. Any baryon asymmetry produced directly in the decay is 
therefore insignificant [12]. However, a nonzero (iV&) will be transferred to baryons 
through the operator (3) if there is an initial particle-antiparticle asymmetry in the 
population of decaying particles. Since e should not be too small, there must exist 
a near maximal asymmetry in the decaying population. 

Such a large asymmetry might appear hard to achieve. However, the coherent 
production of a scalar condensate along a supersymmetric flat direction can give 
rise to a large asymmetry in the condensate [2,6]. Here, flat direction refers 
to a direction in field space on which the perturbative potential vanishes at the 
renormalizable level. Such directions are generic in supersymmetric theories. The 
nonrenormalization theorem protects these directions from being lifted by quantum 
corrections [13]. In the presence of SUSY breaking, a potential can arise though. 
Whether or not a condensate is actually generated along a flat direction depends 
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on the sign of the SUSY breaking soft mass term at early times, m 2 0*0, where 
(f) parameterizes the flat direction. When H m 3 / 2 the finite energy density of 
the universe induces soft parameters along flat directions with a scale set by the 
Hubble constant [6]. If the induced m 2 > 0, the origin is stable, and the large 
expectation values required to form a condensate do not arise. However, if the 
induced m 2 < 0, the origin is unstable and large expectation values can develop. 
In this case, if the flat direction is lifted at order n by a nonrenormalizable operator 
in the superpotential 

W = 13 , n (6) 
then the relevant part of the potential along the flat direction is 

+ m 

where ~ A ~ ^3/2 are s °ft parameters arising from hidden sector SUSY 
breaking, and c ~ a ~ 0(1) are the soft parameters induced by the finite energy 
density [6]. Here the scale M n may in general be (much) different than the scale 
of the operators which allow to decay. For c < the expectation value along 
the flat direction is determined at early times by a balance between the mass 
term and nonrenormalizable terms. If > 0, then when H ~ ^3/2 the origin 
becomes stable and the field begins to oscillate freely with a large initial value. 
However, at just this time since the expectation value of the field is determined 
by a balance between the mass and nonrenormalizable terms, the U(l) violating 
A term necessarily has the same magnitude. Depending on the initial phase of 
the field, the presence of the A term with this magnitude can lead to a near 
maximal asymmetry in the condensate. So if a condensate is produced along a flat 
direction which is lifted by a nonrenormalizable superpotential, it naturally has a 
large asymmetry. For a flat direction made of squark or slepton fields, this is the 
mechanism of baryogenesis proposed by Affleck and Dine [2] . Here however, the 
initial condensate asymmetry is in the (f> field, and is only transferred to baryons 
by decay through the operator (3). 

The final, and perhaps most nontrivial requirement, is that Ulsp + — 1- 
If the decaying particles dominate the energy density at the time of decay, the 
universe is in a matter dominated era at that epoch. However, since ttilsp an d rn§ 
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are the same order this would imply the universe remained matter dominated below 
this temperature. Matter domination from such an early epoch is incompatible 
with nucleosynthesis [7]. While it may have been natural for the condensate to 
dominate the energy density, this is clearly unacceptable [14]. The condensate 
must have a small enough energy density so that matter domination from the relic 
LSPs starts at a temperature of T ~ 5 VLh? eV, where h = H/ (100 km s" 1 Mpc -1 ), 
and H is the Hubble constant now. Assuming critical density, N\£p = 1 in the 
decay, and given the current temperature, the condensate number density at the 
time of decay can be parameterized as 



™LSP 



) 



where s = (27r 2 g*/45)T 3 is the entropy density at the time of decay. This is 
much less than a thermal number density, n/s ~ 1/g*- Now the total density in 
the condensate is determined by the expectation value when the field begins to 
oscillate freely. From (7) the value of the field when oscillations begin (H ~ m 3 ^) 
is 



where a is a constant of order unity. The fractional energy density in the conden- 
sate when oscillations begin is p^/ptot — 0o/^p !■ In an inflationary scenario 
with a reheat temperature low enough to avoid overproducing gravitinos by thermal 
rescatterings, the universe is in an inflaton matter dominated era when H ~ ^3/2 
[6]. So P(j,/ptot stays roughly constant until the inflaton decays. After the inflaton 
decays the condensate density per entropy density is 



^ T R ( m^M 
s ~ m^Ml V /3 



2 

U— 3\ n-2 



(10) 



where Tr is the reheat temperature after inflation. Without additional entropy 
releases n^/s stays constant until the time of decay. So the relic fractional density 
in the condensate is determined by the order at which the flat direction is lifted, and 
the reheat temperature after inflation. For M n ~ M p and n > 6 this is generally 
too large for reasonable reheat temperatures. However for n = 4 and M n / (3 ~ M p 
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, n^/s ~ Tji/Mp. With Tr ~ 10 GeV, this gives just the required condensate 
density. Therefore if the late decaying condensate arises from coherent production 
along a direction which is lifted by a Planck suppressed fourth order term in the 
superpotential, the required relic density naturally arises for reasonable values of 
the reheat temperature after inflation. 

In addition to the required fourth order Planck suppressed term in the super- 
potential, there are in general higher order SUSY breaking terms (in addition to 
the mass term) which are suppressed by the Planck scale. Assuming hidden sector 
SUSY breaking these give a general form for the soft potential of 

V s ^) = ml /2 M 2 p ^/M P ) (11) 

However, just on energetic grounds the nonrenormalizable term in the superpoten- 
tial forces <C M p . The higher order corrections in (11) are therefore unimportant. 
In addition, there are higher order soft terms generated by integrating out fields 
which gain mass at the scale M. These are of the general form 

V s {4>) = ml /2 M 2 g((P/M) (12) 

With M in the range required to give an acceptable decay temperature for 0, these 
higher order terms are less important for H ~ ^3/2 than the terms in (7) with 
n = 4. The higher order SUSY breaking potential terms for therefore do not spoil 
the expectation that the condensate carries a large asymmetry, or the prediction 
for the relic density. 

3. A Toy Model for Baryons and LSPs from Late Decay 

It is easy to build models which satisfy all the requirements outlined in the 
previous section. As an existence proof, consider the following toy model. The flat 
direction required for the coherent production can be parameterized by a singlet 
field 0. In principle this could be a composite field in some sector of the theory, 
but here will be taken to be an elementary singlet for simplicity. The singlet 
should be protected from obtaining a large mass while allowing the operator (3). 
This can be enforced with discrete symmetries. For example, under a Z4 discrete 
R symmetry the superpotential transforms as W — > — W [15]. If and all the 
u and J transform as / — > e" 1 "/ 4 /, where / G (0, u, d), then the operator (3) is 
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allowed while a superpotential mass term, m00, is not allowed. The usual Yukawa 
couplings, X u QH u u, X^QH^d, and XeLH^e, are allowed if the other standard model 
fields transform under the as L — > L, Q — > e Z7r / 4 <5, /i — > e 47F//2 /i, where /i G 
(H u ,Hd,e). The operator (3) must be generated by integrating out particles with 
intermediate scale mass. This can be accomplished in this model by introducing 
a Dirac pair U' and U', with mass, muU'U' , and Yukawa couplings gU'dd, and 
g^ipuU' . These couplings and Dirac mass can be enforced by the transformation 

mass scale m\j could arise from dynamics 
which preserves the discrete symmetry. Integrating out the Dirac pair gives the 
operator (3) with X/M = gg^/mu. The operator (5) is also generated at the scale 
M, but is suppressed by a loop factor compared with (3). So in this model the 
dominant decay mode is — > udd. Finally, the dangerous superpotential coupling 
<f)LH u is restricted by the discrete symmetry. 

The flat direction can be lifted by nonrenormalizable terms in the superpo- 
tential. With the Z4 R symmetry, the lowest order term in the superpotential is 
4 . Such an operator is not generated at the intermediate scale, but presumably 
can arise directly at the Planck scale 

w = wf (13 > 

And, as discussed in the last section, an operator lifted at fourth order in the 
superpotential and suppressed by the Planck scale is precisely what is required to 
give the correct magnitude for the dark matter and baryon densities. In addition 
to the fourth order term in the superpotential, <fi is lifted by a fourth order SUSY 
breaking term in the soft potential generated by integrating out the heavy Dirac 
pair U'U' 

4 2 

V s {<t>) - ^^-(0*0) 2 ^ 

However, as discussed in the previous section, for H ~ ^3/2, terms of this order 
are subdominant compared with (13). 

Acceptable soft SUSY breaking terms can also result in this model. In order to 
allow visible sector gaugino masses the Z4 R symmetry must be broken in the SUSY 
breaking sector to Z2 R parity. For definiteness consider a hidden sector scenario 



9 



in which SUSY breaking is transmitted by Planck suppressed interactions. The 
breaking to R parity can be accomplished with a hidden sector field z which is 
invariant under and breaks SUSY by an auxiliary component expectation value 
(F z ) ~ ^Jmj ) j2Mp. In addition, the soft A term A(3(p 4 /M p , required to generate an 
asymmetry in the <fi condensate, can arise from supergravity interactions, and the 
Kahler potential coupling J d^Ozcj^cj). Dimension 3 soft A terms for standard 
model fields arise from similar couplings. A soft H u Hd scalar mass and weak 
scale /i term can arise from Kahler potential couplings J d 4 6z^ z' H u Hd and 
-jjj- J d 4 9z' H U H^, where z' is a hidden sector field which participates in SUSY 
breaking and transforms as z' — ^ —z' under the discrete symmetry. A weak scale 
mass for the flat direction, m^ 2 0*05 results from supergravity interactions and/or 
Kahler potential couplings with the hidden sector. However, most importantly, 
with the hidden sector couplings sketched above, the Z4 symmetry does not allow 
a soft mass term mg^ 2 00 from Kahler potential couplings, which would violate the 
U(l) carried by <fi. The classical evolution of the condensate at late times therefore 
preserves the asymmetry generated when the coherent oscillations begin. 

So in this model all the requirements are satisfied with a single discrete sym- 
metry. Although the model is perhaps unrealistically simple, it demonstrates that 
all the requirements for baryons and LSPs from late decay of a condensate can be 
met in a technically natural manner. 

4. Conclusions 

In the late decay scenario outlined here, the baryon and LSP densities are 
related by 

- ~ e (15) 

"LSP m LSP 

As discussed above, e < 1 if one LSP results from each decay. Under the assumption 
that the total density is near critical, Q b + ^lsp — 1, a lower limit on the baryon 
density then gives an upper limit on the LSP mass in this scheme. For Q b <C ^lsp> 

m b 

The absolute lower bound on Q b comes from the observed density of luminous 
matter, il b ^ .007. This gives an upper limit of mLSP ^ 140 e GeV. A more 
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stringent upper limit comes from nucleosynthesis. The primordial light element 
abundances depend on the baryon to entropy ratio at the time of nucleosynthesis. 
Comparison of the calculated and observed abundances gives upper and lower limits 
on the baryon density, .01 ^ Vt^h 2 ^ .015 [16]. The nucleosynthesis lower limit 
on the baryon density gives the upper limit ^lsp ^> 100/i 2 e GeV. So the LSP 
is expected to be fairly light in this late decay scenario. For example, in SU(5) 
models for which e < |, with a hubble constant h < .8, the upper limit on the LSP 
mass is mLSP < 30 GeV. 

In addition to the LSPs arising from the late decaying condensate there will 
be a population of LSPs arising from thermal freeze out. However, the low mass 
required for the late decay scenario can give a freeze out density which is well 
below critical. For an LSP which has a sizeable mixture of Higgsino and gaugino 
components, annihilation through s-channel Z exchange is very efficient and leaves 
a very small relic density from freeze out [17,18]. For a mostly Higgsino LSP, 
coannihilation with the other Higgsino states also leads to negligible relic density 
[19]. A small relic density for a light nearly pure gaugino LSP can also result 
from annihilation through t-channel squark and slepton exchange if one of the 
sleptons or squarks are light [18]. So depending on the precise composition of the 
LSP, the late decay can give the dominant contribution to the relic LSP density. 
Independent of the production of a baryon asymmetry, late decay is an interesting 
source of relic LSPs in the low mass regime. In fact, if the LSP was found to be in 
a region of parameter space for which the freeze out density was too small to give 
closure (such as the light Higgsino or mixed Higgsino-gaugino regions) the only 
alternate source for LSPs would be a late decay below the freeze out temperature. 

In conclusion, the ratio of dark matter to baryon density can naturally be 
O(10 — 100) if stable weak scale mass particles and baryons result in roughly equal 
amounts from the late decay of a particle. A natural way in which this can occur 
is for a condensate with a large net asymmetry to decay to R odd combinations of 
standard model fields. Supersymmetric theories with a conserved R parity can in 
principle have all the ingredients to realize this scenario. 

I would like to thank R. Brandenberger, M. Dine, J. Primack, L. Randall, and 
U. Sarid for useful discussions, and J. Frieman, R. Malaney, and G. Steigman for 
discussions about nucleosynthesis. 
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